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Abstract 



By calculating the contribution of the it — it three-body force to the three- 
nucleon binding energy in terms of the ttN amplitude using perturbation 
theory, we are able to determine the importance of the energy dependence 
and the contribution of the different partial waves of the ttN amplitude to the 
three-nucleon force. A separable representation of the non-pole ttN amplitude 
allows us to write the three-nucleon force in terms of the amplitude for A^A^ — > 
NN*, propagation of the NNN* system, and the amplitude for A^A^* 
NN, with A^* being the vrA^ quasi-particle amplitude in a given state. The 
division of the ttN amplitude into a pole and non-pole gives a procedure 
for the determination of the ttNN form factor within the model. The total 
contribution of the three-body force to the binding energy of the triton for 
the separable approximation to the Paris nucleon-nucleon potential (PEST) 
is found to be very small mainly as a result of the energy dependence of the 
ttN amplitude, the cancellation between the S- and P-wave ttN amplitudes, 
and the soft ttNN form factor. 
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I. INTRODUCTION 



The discrepancy between the results of the exact calculations for the binding energy 
of the triton using a number of realistic nucleon-nucleon potentials and the experimental 
value of 8.48 MeV, has been an outstanding problem in nuclear physics for a number of 
years |jl|]. A commonly accepted solution has been the introduction of a three-nucleon force 
that will bridge the gap between the calculated binding energy 0-^, based on a two-body 
interaction, and the experimental binding energy. The origin of such a three-body force is 
partly the result of the fact that the nucleons are treated as point particles interacting via 
a two-body meson exchange potential that is often assumed to be local. The fact that the 
formal division of the interaction into a two- and three-body part is not unique, given the 
on-shell two-body data ||^, suggests that the contribution of the three-body force is partly 
determined by the definition of the two-nucleon interaction. Thus under ideal conditions 
the division of the interaction, in the three-nucleon system, between a two- and three-body 
force will require a consistent formulation of these two potentials within a meson-nucleon 
theory. 

In the absence of such a formulation one may assume that a meson-nucleon theory should 
give the correct binding energy for the three-nucleon system, in which case the three-body 
force is by definition that force which when added to the chosen two-nucleon force will give 
the three-nucleon binding energy l^HH- ^ second approach is to assume that the three- 
nucleon force is the result of meson exchanges that are possible only when the number of 
nucleons is greater than two. In this second approach one expects the dominant mechanism 
to be one in which one nucleon emits a meson that scatters off a second nucleon and then gets 
absorbed on the third nucleon, see Fig. |l|. In this case the three-nucleon force is determined 
by the off-shell meson-nucleon amplitude that goes into the calculation of the diagram in 
Fig. |1|. In the present investigation we will consider the second approach involving a vr 
meson exchange. In particular, we will examine the role of the energy dependence of this 
vriV amplitude on the contribution of this three-body force to the binding energy of the 
triton. We will also examine the relative contribution of the different vriV partial waves to 
this three-body force. 

Over the past ten years three approaches have been developed to determine a three- 
nucleon interaction from ttA^ dynamics, (i) The Tucson-Melbourne (TM) [|Tn],|n| three- 
nucleon potential is based on the idea that the off-mass shell vrA^ amplitude should satisfy 
current algebra constraints and the soft pion theorems [0. These constraints allow a 
covariant parameterization of the off-mass shell amplitude. To be consistent with the meson 
exchange A^A^ interaction, the vrA^ amplitude is expanded in powers of where q is the 
pion momentum and rriN the nucleon mass. This gives a vrA^ amplitude that includes both S- 
and P-wave scattering, but where the energy dependence is reduced to z/ = (s — u)/Amjq = 
{q' + q) ■ {p' + p)/4m]si = 0, as a result of the expansion in Here, s and u are the 

usual Mandelstam variables while q (g') and p {p') are the initial (final) four momentum 
of the pion and nucleon respectively. The ttNN form factor is constructed to satisfy the 
Goldberger-Treiman relation p!3|JT^ . Although the original TM potential included only the 
vr — vr three-body force, ir — p and p — p contributions have recently been included JlS 



(ii) A similar approach is to assume that the vrA^ dynamics is determined by an effective 
chiral Lagrangian |[I6H20| which when used to calculate the nN amplitude at the tree level, 
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will give an effective three-nucleon force given by the diagram in Fig. |l]. The evaluation of 
the ttN amplitude at the tree level gives rise to an energy independent vrA^ amplitude, and 
therefore a three-body force. Both approaches (i) and (ii) give similar results, and emphasize 
the chiral symmetry of the vrA^ amplitude, (iii) The iV — A coupled channel approach [p]-|25[| 
takes as its starting point the fact that the vrA^ amplitude is dominated at medium energies 
by the A(1230) or the P33 partial wave. This dominance of the A suggests that we could 
extend our Hilbert space to include not only the nucleon, but also the A as an excited 
state of the nucleon. The approach of treating the N and A on equal footing effectively 
includes, in a consistent manner, that part of the original three-body force corresponding to 
T^N scattering in the P33 channel, or at least the resonance part of it [Q. The advantage 
of this approach is that now we can construct the two-body and the three-body forces 
with some consistency, to the extent that the A^A^ — A^A transition potential used in the 
two-body interaction can also be used to generate the three-body potential. The inclusion 
of coupling between the A^A^, A^A and AA channels allows a consistent treatment of the 



BB and BBB system, where B = N,A [25|. In this approach, since the vrA^ amplitude 



is basically approximated by vrA^ ^ A ^ ttN, the energy dependence of the total vrA^ 
amplitude is completely determined by the energy dependence of the P33 channel which 
has the A resonance. Although the A^A^ and transition potentials can be local and energy 
independent, the effective three-body force in this model is energy dependent and this energy 
dependence is determined by the A resonance, i.e. the P33 amplitude. More recently there 
have been extensions of the A^ — A coupled channel approach that have included the S'-wave 
component of the TM potential . Also Pena et al. have examined the importance of 
the coupling of the A to the vrA^ channel. This latter calculation gives a vrA^ amplitude in 
the P33 channel that fits the phase shifts and has an energy dependent mass and width for 
the A. The inclusion of this energy dependence in the A mass and width does not effect 
the final result appreciably. 

The questions which arise from the above approaches to the vrA^ dynamics that go into 
the derivation of the three-nucleon force and its contribution to the binding energy of the 
three-nucleon system are: (i) What are the contributions of the different vrA^ partial waves to 
the three-body force? (ii) If there is any cancellation between the contribution of the different 
partial waves, should the P33 partial wave be treated via the N — A coupled channel while the 
rest of the vrA^ amplitude gives rise to a three-body force as depicted in Fig. ^ (iii) Would 
a cancellation between the different vrA^ channels be sensitive to the energy dependence of 
the amplitudes? The main aim of this investigation is to examine these questions. 

To motivate our interest in the importance of the energy dependence of the vrA^ am- 
plitude, let us examine the role of the A^A^ amplitude in calculating the binding energy 
of the three-nucleon system within the framework of the Faddeev equations. Here we ob- 
serve that we require the fully off-energy- shell A^A^ amplitude in a given partial wave a, 
t^'^ {k, k'; Enn), for all energies in the range —00 < -Eatat < —Et, where Et is the three- 
nucleon binding energy, see Fig. ^. The fact that we need this amplitude over the full 
specified energy domain is a result of the fact that in the three-nucleon system, the total 
energy is fixed at E = —Et, and the spectator particle can have any energy from zero to 
00. To consider the contribution of Fig. |1] and in particular the energy dependence of the 
ttN amplitude, we consider the NNN — ttNNN equations [p8| , p9| which are an extension of 
the A^A^ — T^NN equations to the A = 3 sector. Within the framework of the above ttNN 
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dynamics, to calculate the contribution of the three-body force as defined in Fig. 1 to the 
three-nucleon binding, we must determine to know the fully off-energy-shell ttN amplitudes 
30| , f^^i^k, k'; -EttAt), for all energies in the range —oo < E^^^ < {uin — Ex), where in this 



case we have included the rest mass of the nucleon and pion in the tiN energy E^^^- In other 
words we have to calculate the fully off-energy- shell vrA^ amplitude for all energies from —oo 
to Et below the nucleon pole ( see Fig. |^). Thus to calculate a three-body force which is 
defined in terms of the ttN amplitude, we need to know this amplitude off-shell for energies 
far below {{m.^ + Et) ~ 145 MeV), the nN threshold. Here the need for the ttN amplitude 
over this energy domain is a result of the fact that the total energy is still —E^, but now we 
have two spectator nucleons with kinetic energy between zero and oo. At these energies it is 
not clear that the A(1230) is dominant or that the P-wave amplitudes are more important 
than the S'-wave amplitudes. Clearly, the threshold behavior of the amplitude is essential, 
but we need to know the nN partial wave amplitudes over a wide range of energies above 
threshold, if we are to extrapolate these amplitudes to the energies required in calculating 
the three-body force. This suggests that we have to fit the energy dependence of the ex- 
perimental vrA^ phase shifts in order to improve the accuracy of the extrapolation in the 
energy. 

In Sec. |I| we will derive the three-body force using the coupled channel method and we 
will discuss approximation which we take to avoid the dressing problem In Sec. ^ 



we will discuss a parameterization for the vrA^ amplitude, and in particular will discuss the 
division of this amplitude into a pole and non-pole component. It is only the non-pole part 
of the ttN amplitude that goes into the nN scattering that generates the three-body force. 
This division is essential if we are to avoid double counting. We will choose the irNN form 
factor as the residue of the pole term in the vrA^ amplitude to keep the consistency with the 
nN formulation. We will then proceed in Sec. |rV| to a discussion of the vr — vr three-body 
force given in Fig. ^ in terms of our parameterization of the vrA^ amplitude. As a result of 
using a separable representation for the vrA^ amplitude, we find that the three-body force is 
reduced to the product of the amplitude for A^A^ — > NN*, followed by the propagation of 
the NNN* system, and the final transition amplitude for A^A"* — »• A^A^, where a runs over 
all 7cN partial waves. This tt — vr three-body force is employed in Sec. |V] to calculate, in the 
Born approximation, the contribution of this force to the binding energy of the three-nucleon 
system. The three-nucleon wave function will be calculated by using Paris (PEST) potential 
|3^ . Since the Paris potential does not have any energy dependence, our approach for the 
three-nucleon interaction is inconsistent with the two- nucleon interaction. To that extent, 
the overall magnitude of the resultant three-body force contribution may not have great 
significance. However, we would like to examine the relative contribution of the different 
vrA^ partial waves, and see how this contribution is sensitive to the energy dependence of 
the irN amplitude in the subthreshold region. In particular, we would like to demonstrate 
how the energy dependence of both the non-pole amplitude and the nNN form factor play 
a role in the importance of this three-nucleon force. Here we will find that the energy 
dependence of the non-pole amplitude has a significant effect on the contribution of the 
three-body force. Furthermore, we will demonstrate that there is a cancellation between 
the contribution of the different vrA^ partial waves. This cancellation turns out to be very 
sensitive to the approximations used. Finally, in Sec. we will conclude our discussion by 
considering some open questions that can influence our final results. 
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II. COUPLED CHANNEL FORMULATION OF THREE-BODY FORCE 



In this section, we will establish the approximation involved in writing Fig. |I| as the 
lowest order contribution to the binding energy of the three-nucleon system from a tt-tt 



three-body force. Our starting point is the Hamiltonian of Mizutani and Koltun [p3|l 



n = lC + VNN + V^N + A + A^ , (2.1) 

where /C is the kinetic energy of the nucleons and pions, while Vnn and Vt^n are the A^A^ 
interaction in the absence of one pion exchange, and the ttN interaction in the absence of 
the s-channel nucleon pole, respectively. In Eq. ( p.l| ), A is the pion absorption vertex while 
A^ is the corresponding production vertex. The Schrodinger equation for this Hamiltonian 
is 

H\^)=E\m). (2.2) 

Clearly, the operators A and A^ change the number of pions. To that extent, the present 
Hamiltonian has some of the features of a field theory in that the number of pions is not fixed 
while the number of nucleons is fixed. As a first approximation, we restrict the Hilbert space 
to [nN) and 7i{nN) systems only, where n is the number of nucleons, which is conserved. 
We define the Feshbach [^] projection operators onto [nN] and n{nN) spaces as P and Q, 



respectively. By using these projection operators and by assuming that this truncated space 
is complete, i.e. P + Q = 1, the wave functions of the (nN) and T[{nN) components ( P^ 
and , respectively) are solutions of the equations 



E-Hpp- HpQ Hqp) \Pm) = , (2.3) 

-C/ — linn ' 



and 



E-Eqq- Hqp Hpq) m) = . (2.4) 

11/ — npp 



In writing Eqs. (|2.3|) and (|2.4|) we have made use of the definitions. 



Hpp = PUP , HpQ = PHQ , 

Hqp = QHP , and Hqq = QHQ . (2.5) 

The Green's function for the 7T{nN) system, {E — Hqq)^^, which including the vrA^ and A^A^ 
interactions, can be written as 

^ GlAE) + G^XmiE) T{E) G^^JE) , (2.6) 



£j _ ff^Q ^ T7(nN)\ J ' ^ iT{nN)\ J \ I ^ TT{nNy 

where G^Xj^-^{E) is the free 7r{nN) propagator. The corresponding T- matrix, T{E), can be 
written as 

T{E) = UE) + UE) Gi%(E) U^,{E) Gi%(E) t,{E) , (2.7) 

a ab 
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where Uab{E) is the AGS amphtude [^] for the (n + 1) particle system and satisfies a set 
of coupled equations. Here, ta{E) is the vrA^ or A^A^ amplitudes resulting from the potential 
Vnn or Vttat, respectively. 



We will first consider the simplest case of n = 1 [^. Here, on the one hand, Eq.( p.3|) 
can be regarded as the equation for the dressed nucleon with the self-energy, defined 
to be 



= HpQ 



1 



E-Hr 



QP ■ 



(2.8) 



Therefore the dressed nucleon propagator, Gj\f{E), can be written as 

\ -1 



Gn{E) 



E — Hpp — H 



1 



E-Hr 



QP 



(2.9) 



with Hpp, the bare mass of the nucleon. On the other hand, Eq. ( |2.4| ) describes vrA^ 
scattering. Now since Hqq = Kqq + VttN, with Kqq the kinetic energy of the vrA^ system, the 
irN interaction in Eq. ( ^.4[ ) is the sum of two contributions. The first is the irN interaction 
in the original Hamiltonian, V^tn, while the second term, Hqp {E — Hpp)~^ HpQ, results 
from the coupling of the ttN channel to the A^ channel. The ttN T- matrix, T^Ar(i?), in this 
case is the sum of two terms, and is given by 



TME) = Fqp{E) Gn{E) Fpq{E) + tB{E) 



where 



Epq{E) = HpQ 



1 + 



E-Kqc 

and tB^E) is the solution to the two-body equation 

tB{E) = VQQ + VQQ ^ 



■tB{E) 



(2.10) 



(2.11) 



E-Kc 



tB{E) 



(2.12) 



with Vqq = VnN- Here we note that the first term on the right hand side of Eq. ( |2.1U| ) 
corresponds to the process vrA^ —^N^ nN, and as a result the pion absorption (emission) 
vertex is given by Fpq{E) {Fqp{E)). To establish that this term has the nucleon pole 
contribution to the vrA^ amplitude, we note that the Green's function, Gm{E), given in 
Eq. ( p.9|) , has a pole at the nucleon mass, and therefore can be written as 



Gn{E) 



1 



Z2{E) 

E - Hpp - S(E) E-tun 



(2.13) 



where Z2{E) is the wave function renormalization. Making use of this result in Eq. ( p.lO| ) 
allows us to write the ttN T-matrix as 



TME) = f^LiE) 



&N{E)+tB{E) 



E — rriM 

where the renormalized nNN form factor, which is energy dependent, is given by 



(2.14) 
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f^^^iE) = Zi{E)FpQiE) . (2.15) 

In this formulation the normahzed physical nucleon wave function, which is a solution to 
Eq. ( p.3|) is given by 

^ N = zl (niN) P"^ ■ (2.16) 

We consider next the case of n > 1. Since we have restricted the Hilbert space to include 
n nucleons, and up to one pion only, this truncation effects the dressing of the nucleons for the 
case n > 1. From Fig. ^, we observe that the nucleons can only be dressed separately after 
absorbing an initial pion and before emitting the final pion, because of the limitation imposed 
on the Hilbert space. We can neither include the nucleon dressing before the pion absorption 
nor after the final pion emission. This incomplete dressing makes the renormalized nNN 
form factor, /^^(i?), smaller To overcome this problem we need to guarantee that 



all the nucleons are fully dressed at the same time, and this dressing is on both sides of 
the TrNN-Yertex ^7^, see Fig. |^. We will avoid this problem in the present investigation in 
the following manner. In the NN — ttNN equations the nucleon dressing was introduced 
to satisfy the unitarity of the A^A^ amplitude above the pion production threshold. In the 
present investigation however, we are considering a bound state problem which is below 
the threshold for pion production. We therefore expect that the nucleon dressing will not 
be essential, and does not have to be included explicitly to satisfy unitarity. We therefore 
assume that each nucleon line has been renormalized to give the physical nucleon mass and 
correct ttNN coupling constant, i.e.; (i) The nucleon propagator is given by {E — mjv)~^. 
(ii) The ttNN vertex is given by the renormalized vertex function, /^^(i?). 

Using the above assumptions, the interaction term resulting from the coupling of the 
NNN and irNNN Hilbert space in Eq. ( p.3|) can be rewritten as 

HpQ 

where the series is generated by iterating the equation for Uab in Eq. (|2.7|) , and making use of 
the resultant nNNN T-matrix in Eq. (|2.6| ). In particular, if we go to third order in the vrA^ 
amplitude ta, the vrA^A^ form factors in V3B get dressed as detailed in Eqs. (|2.14|) and ( p.l5[) . 



^ ^rr Hqp = VoPE + Vd^sp + V3B + " " " , (2.17) 



Each of the terms in Eq. ( p.l7| ) is shown schematically in Fig. ^ We note that Vqpe and 
Vdisp are two-body operators. Since we will be calculating perturbatively the contribution 
of the three-body force to the binding energy of '^H, the three nucleon wave function used 
results from the solution of the Faddeev equations for a given nucleon-nucleon interaction. 
Assuming that the nucleon-nucleon interaction includes Vqpe and Vdisp, we omit these terms 
from Eq. (|2.17|) . This leaves us with the three- nucleon force Vsb, illustrated in Figs, ^(c) 



and 0, which gives the lowest order correction to the binding energy of the three-nucleon 
system. The detail of the numerical framework of the vrA^ scattering will be presented in 
next section. 



III. THE ttN amplitude 

To include the full energy dependence of the vrA^ amplitude into the determination of 
the 7T — 7c three-nucleon interaction and its contribution to the binding energy of the triton. 
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we have to: (i) Remove the nucleon pole contribution from the vrA^ amplitude to avoid 
double counting in the three nucleon force calculation, (ii) Define a vriViV form factor for the 
emission and absorption vertices in Fig. |1] that is consistent with the vrA^ amplitude used and 
the scattering data, (iii) Treat the nucleon in the vrA^ system using non-relativistic kinematics 
to maintain consistency between the NNN and vrA^ systems since the nucleons in the three- 
nucleon system are treated non-relativistically. From the previous section we observe the 
first two conditions can be satisfied if we choose a formulation of the vrA^ scattering problem 
that is motivated by a Hamiltonian that includes a ttNN vertex and a ttN interaction, e.g. 

m 



the Cloudy Bag Model 



Thus a choice for the vrA^ potential, motivated by the lowest 



order contribution to the amplitude and based on a Lagrangian of the form suggested by the 



Cloudy Bag Model with volume coupling ||3^, consists of an s-channel nucleon pole diagram, 
Fig. 1^, and the cross diagram and contact term. Fig. ^. 

For the present investigation, to simplify the parameterization of the ttN amplitude, we 
replace the cross diagram and contact term in Fig. |^ by a one term separable potential in 
each partial wave. This allows us to write the vrA^ potential in a given partial wave a as 



E — niQ 



h{k')+g^{k) Kg^{k') 



(3.1) 



where the first term corresponds to the nucleon pole diagram with a bare nucleon mass of 
mo, and a bare form factor fo{k). Since we are using non-relativistic kinematics for the 
nucleon, this s-channel pole diagram contributes to the Pn partial wave only. As a result, 
the potential in all partial waves other than the Pu channel is given by the second term on 
the right hand side of Eq. (p.l|). This separable vrA^ potential gives an amplitude of the form 



taik, k'; E) = g^{k) Ta{E) g^{k') 



for a ^ P] 



11 



where 



with 



T^{E) = - {ga\G^N{E)\g^) 



-1 



OO 

{ga\G.nN{.E)\ga) = j dkk" 



E — Uk — r 



k' 



(3.2) 



(3.3) 



(3.4) 



- rriN 



This amplitude is a solution of the Lippmann-Schwinger equation with the pion energy 
treated relativistic since Uk - 



'k"^ + m^. The strength of the potential A^, and the form 
factor ga{k) are adjusted to fit the experimental vrA^ phase shifts in all S- and P- waves other 
than the Pn channel. For the present investigation we will use the parameterization used 
by Thomas for -rid scattering. The form factor used for Sn and 5*31 vrA^ partial waves is: 



ga{k) 



Si 



k^ + al 



+ 



and 



9a{k) 



Si 



k^ + aV 



s,.e 



k^ + ai (F + ai)2 



(3.5) 



(3.6) 
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for Pi3, P31 and P33. This parameterization has also been used extensively in the NN — ttNN 



calculation for vrrf scattering and Tid pp reactions . 

The Pii channel plays the important role in this analysis of the three-nucleon interaction 
as it has the nucleon pole contribution that needs to be removed to avoid double counting. 
It also has the information about the tcNN form factor which is defined as the residue of the 
off-shell nN amplitude at the nucleon pole. In this way we can subtract the nucleon pole 
contribution to the irN amplitude, and extract a ttNN form factor while maintaining a fit 
to the TcN scattering data in this channel. Since the potential in this channel is the sum of 
two contributions, see Eq. (^.11), we can write the corresponding amplitude as a solution of 
the Lippmann-Schwinger equation using a two-potential theory to be |43] 



Uk, k'- E) = f{k; E) ^ f{k'- E) + g^{k) t^{E) g^{k') , (3.7) 

where the dressed tcNN form factor, f{k; E), is given by 

f{k;E) = fo{k)+g^{k)T^{E){g^\GME)\fo) , (3.8) 

and a in Eq. ( |3.7| ) and ( |3.8| ) refers to the Pu channel. In Eq. ( ^.71 ), the mass renormalization 
factor S(i?) is given by 

J:{E) = {fo\GME)\f{E)) 

= {fo\G^N{E)\fo) + {fo\GME)\g) r{E) {g\GME)\h) , (3.9) 

where we have dropped the channel label a with the understanding that we are considering 
the Pu partial wave only. 

In writing the amplitude in the Pu channel ( Eq. ( |3.7| ) ) as the sum of a pole and 
a background term, we are able to define that part of the nN amplitude which will be 
included in the evaluation of the three-body force. At the same time, we can determine the 
ttNN form factor that is required for the pion emission and absorption vertices in Fig. |l]. To 
establish that this nNN form factor gives the correct ttNN coupling constant as the residue 
of the vrA^ amplitude at the nucleon pole, we expand S(-E) about the physical nucleon mass 
as 



S(E) = S(mjv) + {E- mjv) Si(m^) + {E - m^)^ j:^{E) . (3.10) 
If we now fix the bare mass tuq such that 

mo S(mAf) = m-AT , (3.11) 
we can write the pole amplitude as 

m E) - — ^— - f{k'- E) = nk- E) — ^ nk'- e) (3.12) 

E — mo - 1^{E) E — niN 

where the renormalized tcNN form factor f^{k] E) is defined as 
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and the wave function renormalization constant Z2 is given by 



Z2 = [1 - Si(mAr)] 



-1 



with 



Sf(E) = Z2S,(E) for « = 1,2 



(3.14) 



(3.15) 



In this way we have defined the renormahzed irNN form factor, f^{k; E), which will be used 
for the emission and absorption vertices in Fig. More important is the fact that this form 
factor is constrained by the vrA^ phase shifts in the Pu channel and the requirement that we 
have the correct irNN coupling constant. We note at this point that the renormalized nNN 
form factor f^{k; E) is energy dependent and that this energy dependence is determined 
by unitarity through Ta{E), and has to be included if we are to fit the phase shifts in this 
channel. 

To determine the nNN coupling constant resulting from the above formulation of vrA^ 
scattering, we need to compare our results for the pole amplitude with the corresponding 
Feynman diagram for the Lagrangian with the pseudoscalar coupling, i.e. 



(3.16) 



where the coupling constant 5^0 is made a function of the momentum. In Eq. (|3.16|) , r is the 
Pauli isospin matrix, (f) is the pion field, and u is the usual Dirac spinor. This interaction 
Lagrangian allows us to determine the invariant amplitude corresponding to the s-channel 
nucleon pole diagram [^, and the corresponding S'-matrix. Making use of the relation 



between the S'-matrix and the T-matrix [46|, which is a solution of the Lippmann-Schwinger 
equation, we can calculate the ttNN coupling constant as a result of the relation between 
the T-matrix, t{k, k'; E), and the invariant amplitude T{k, k'; E) as 



t{k, k'- E) = C{k; E) T{k, k'; E) C{k'; E) 



where 



C(k-E) 



\ {27r)^uJkekiE + ujk + £k] 



(3.17) 



(3.18) 



with Ek = yk'^ + m%. The coupling constant is now defined as the residue of the invariant 
irN amplitude at the nucleon pole with all the legs of the nNN vertex on-mass-shell. This 
corresponds to taking E = mj^ and k = ko, where 

l--4l . (3.19) 



7 2 2 



This definition allows us to write the nNN coupling constant f-rrNN in terms of the renor- 
malized TcNN form factor f^{k; E) as 



9l{k) 



ml mNiek + rriN) 



Am% 67i[C{k,mN)f 



f^{k;mN) 
k 



(3.20) 
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where 



^/'iV7v(M = 0.079. (3.21) 



For the present investigation we make use of the Pu nN parameterization of McLeod and 
Afnan where the bare ttNN form factor fo{k) is taken to be 



^o(^) = ^ (1.2 ^ 2,.o ' (3.22) 



(3.23) 



while for the background separable potential form factor, we take 



Cl ^ C2{kY' 



k^ + Pf {k^ + (3lf 

This choice of form factor is basically the same as that used by Thomas ( Eqs.( p.5|) and 
( p.6|) ). The factor of ^JuJ^ was introduced to get the covariant phase space to determine 
the coupling constant at the nucleon pole with all legs on-mass-shell. In Table | we present 
two parameterizations of the Pu amplitude corresponding to a monopole { = 1 ) 
or a dipole ( ng = 2 ) bare ttNN form factor. The parameters were adjusted to give the 
phase shifts below the pion production threshold, the position of the nucleon pole, and 
the ttNN coupling constant of 0.079. In Table |I| we present the scattering volume an in 
this channel, the wave function renormalization Z2 and the value of the coupling constant 
at A; = 0. We note that the values in Tables | and || are the corrected values for the 
parameters of these potentials and their predictions. We note that the renormalized form 
factor f^{k;E) is substantially different from the bare form factor fo{k) due to the wave 
function renormalization Z2 and the contribution of the non-pole amplitude to the form 
factor dressing (see Eq.( |3.8| ) and ( |3.13| )). Finally, the value of the coupling constant at 
k = 0, when compared with the value at k = ko, can be used as a measure of the deviation 
from the Goldberger-Treiman relation. 

The present definition of the irNN form factor is different from that commonly used 
in NN potentials, and three-nucleon force calculations. However it is consistent with the 
formulation of vrA^ scattering where the ttN amplitude is a solution to a two-body equation 
47| . Traditionally the ttNN form factor, introduced as a cutoff in the NN amplitude, is a 



function of the pion momentum only. This is a result of taking both nucleons in the vertex 
on-mass-shell. However in a non-relativistic or time ordered theory, intermediate particles 
are off the energy shell. As a result, the irNN form factor becomes a function of the energy 
and the relative momentum. The energy dependence of the nNN form factor is the result of 
the dressing. This dressing is necessitated by the requirement that the full nN amplitude, 
even in the Pn channel, should be a solution of the Lippmann-Schwinger equation. In 
this way we can maintain consistency with both of the treatments for the nN and NN 
amplitudes as solutions of two-body equations when used in the three-nucleon system. 

IV. THE THREE-BODY FORCE 

Having defined our vrA^ amplitude, and in particular how the nucleon pole is subtracted 
from the vrA^ amplitude to give us a nNN form factor that is both momentum and energy 
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dependent, we turn our attention to our definition of tlie tliree-body force as given in Fig. |T|. 
At this stage it is important to note that this definition of the three-body force does not 
include all possible pion exchange diagrams that are not included in the nucleon-nucleon 
interaction. However, we expect the diagram in Fig. |I] to give the dominant contribution to 
the three-body force. 

To evaluate the diagram in Fig. |l| we introduce Jacobi variables in the irNNN center of 
mass. These are defined in Fig. ^. This choice for the momenta will allow us to take matrix 
elements of the three-body force between three-nucleon wave functions resulting from the 
solution of the Faddeev equation in momentum space for a given two-nucleon interaction. 
At the same time we will be able to include the vrA^ amplitudes defined in Sec. ^TT|with their 
full energy dependence with no approximations. The momenta in the initial state in Fig. ^ 
are 



qa = -ks (4.1) 

mAr(k^ + k'J - (m7v + m^)k2 
P3 = 77^ . ^ (4.2) 

where k'^^, k2, ka and k^ are the momenta of the three nucleon and pion after the pion 
emission vertex. Here, Q3 is the tcN relative momentum for the pion production vertex, 
while Q3 is the relative vrA^ momentum in the vrA^ amplitude. In a similar manner we 
can define all the Jacobi momenta before and after the nN scattering and before the pion 
absorption. In this way all momenta are defined in terms of the initial and final momenta of 
the three nucleons. At this stage we should point out that in a non-relativistic theory, the 
ttNN vertex is not Galilean invariant since mass is not conserved. As a result, the relative 
momenta P3 and qa are not the same before and after the pion emission. For practical 
calculations we will assume that not to be the case. In other words the relative momenta 
Ps and qs are those used in the three-nucleon wave function resulting from the solution of 
the Faddeev equations. 

The pion absorption and emission vertices in Fig. ^ have the form factor f^{Qi; Ei), i = 
1, 3 respectively, where the energy Ei is the energy available to the vrA^ system, and is given 
by 

E, = E + mr,-^-^ , (4.4) 

where E = —Et is the total energy of the system not including rest masses, and the reduced 
masses /xi and /X2 are defined by the relations 

11 1 J 1 1 1 /. 

— = + 7 , and — = + . (4.5) 

In a similar manner, the vriV scattering in channel a in Fig. ^ is represented by the amplitude 
ta{Qi, Q3', E2) which is given by 

t„(Q;, Q'3; E2) = 9a{Q[)Ta{E2)9a{Q3) , (4.6) 
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where the energy available to the nN system, i?2, is given by 

E, = E + m^-^-^ = E + m^-^-^ . (4.7) 
2/i2 2/ii 2/^2 2/ii 

In this way we have made use of the general structure of the nN amplitude in terms of a 
one-particle reducible (the s-channel pole amplitude) and the one-particle irreducible (the 
non-pole amplitude) to determine the three-body force contribution from Fig. |l]. Although 
we have used a separable potential for the non-pole amplitude, there is no reason why we 
could not have made use of the non-pole contribution from a chiral Lagrangian or a 



separable approximation to such a chiral nN amplitude other than the fact that this 
would have imposed an additional complexity to the evaluation of such an amplitude. 

Having defined the basic ingredients required to calculate the contribution from the 
diagram in Fig. |l|, we turn our attention to the practical problem of calculating the overall 
contribution of the diagram in Fig. |T]. The contribution from the process whereby nucleon 
1 emits a pion that scatters off nucleon 2 in the vrA^ channel a, and then gets absorbed on 
nucleon 3 is given by the expression 

X {gaW{p'3q3\G^NNNmP3q3){Q3\f''m) , (4.8) 

where /?i and (3^ represent the set of quantum numbers that label the three-body channels 
in the final and initial states, respectively. The coefficient F^^^^ is a factor determined by 
the transformation between the different Jacobi momenta in the three-nucleon system. The 
four-body Green's function Gt^nnnIE) can be written in terms of the nNN Green's function, 

G-rrNN, aS 



GTtNNN{E) — GttNN ( E — — 




2 



= (E,-m^-^-sjQl + mlj , (2 = 1,3). (4.9) 

This allows us to employ the methods developed for pion exchange in the NN — nNN 
problem and to write the partial wave projection of the process whereby nucleon 1 emits a 
pion that will scatter off nucleon 2 in channel a as 

Zfk (p's,P3;E-^] = {9am{p'3q3\G^NNN{E)\psq3){Q3\f''{E^)) 



3 '^3 



2/^2, 



AL,a,b e^-a+b I 1+a-b a+b 
2^ Ps Ps P2 

L,a,b 



where ij^ is the relative vriV orbital angular momentum in the vrA^ amplitude. Here /^g 
gives the set of quantum numbers for the coupling scheme [(7rA''i)Ar, N2] resulting in a total 
angular momentum jg and total isospin t'^. In a similar manner a'^ = {a, 73) gives the set of 
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quantum numbers for the couphng scheme [{'n'N2)a, Ni] giving rise to the same total angular 
momentum jg and total isospin tg. The quantum numbers of Ni are given by 73 = (£3, S'^), 
where £3 is the orbital angular momentum of Ni relative to the center of mass of {itN2), 
and 5*3 ( the corresponding channel spin) is the sum of the total angular momentum (j^r) of 
(7rA^2), and the spin (si) of nucleon Ni, i.e., S3 = + Si. In writing the above expression 
for the one pion exchange amplitude, we have made use of the fact that both Q3 and Q3 
can be written in terms of the momenta ps and P3 as 

Q3 = P3 + P2P3, Q3 = P3 + P2P3, (4.11) 

where 

P. = . (4.12) 

rriN + m^n- 

and y = Ps ■ P3. The coefficients A^f^, are those used in the partial wave expansion of 
the Faddeev equation for a separable potential |^^. For the case when the vrA^ channel a 
corresponds to the P33 partial wave, Eq. ( |4.1CI| ) gives the (^'3,^3) partial wave projection of 
the NN — NA transition potential. In general, since we have restricted our analysis to the 
case of separable non-pole vrA^ amplitudes, we can interpret the irN amplitude in each partial 
wave to be dominated by an N* in which case Z V 'i, can be considered as the transition 
potential for A^A^ NN*. The difference between the traditional A^A^ — A^A transition 
potential and the above result in Eq. ( |4.1CI| ) is the fact that the present transition potential 



is energy dependent and the parameters of the potential are determined by the vrA^ data 
rather than by the A^A^ data. 

This interpretation of Z^J'^i^, as a partial wave projected transition potential will allow 
us to regards the diagram in Fig. |l] to correspond to an initial state of three nucleons with 
nucleons 1 and 2 in the channel (3'^ going to two nucleons plus an A^*, followed by the A^* 
coupling to nucleon 3 to form a final state of nucleon 2 and 3 in channel (3[ with nucleon 1 
as spectator. This corresponds to the exchange of an A^* and will allow us to partial wave 
expand the three-body force in a similar manner to the expansion of the one pion exchange 
diagram treated above. The resultant partial wave expansion of the diagram in Fig. ^ is 
given as 

r^{E2)Z'^^l,, {^^,V^-E-^ , (4.13) 

where the phase, R = Jtt + tvr + ^3 + 5*3 + tg, results from changing the coupling scheme 
to maintain consistency with the definition of the Jacobi coordinates. Here, {jnytn) are the 
total angular momentum and isospin of the pion with nucleon 2, i.e., they define the vrA^ 
channel a, while £-(2 = 1,3) are the orbital angular momenta corresponding to the momenta 
p-, e.g. £3 is the orbital angular momentum of A^i relative to the center of mass of (7rA'2). 
In this case, we have written the momenta pg and p'^^ in terms of the momenta q3 and qi as 
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Ps = qi + pi qs , p'l = qs + pi qi , 



(4.14) 



where 

and X = qi ■ qs. Here again the coefficient ^^{"/^g is the coefficient resulting from the angular 
momentum recoupling, i.e., [(A^^iVa)/?^, ^"3] ^ 'Ni + N; + N3 ^ [(iV2iV3)/3i , iVi]- 

To compare the three-body force given in Eq. ([4.13|) with that resulting from the 
A^A^ — A^A coupled channel approach [^], we note that the expression in Eq. ( [4.13|) ba- 



sically consists of a transition potential A^A^ NN* followed by the propagation of the A^* 
quasiparticle and finally the transition potential A^A^* —>■ NN. Making use of the definition 
of Ta{E2), ( Eqs. (|0|) and p.4|) ), we can write this quasi-particle propagator for the case 



when a refers to the P33 channel as 

r,(E2) = [(^2 - E^){g^\GME2)G^N{EA)\ga)Y 

[{ga\G-„N{E2)G.„N{EA)\ga)] ^ 



I3 _ P'l 



E + - - 2^ - 2^, 



(4.16) 



In writing this equation we have made use of the fact that t^^{Eis) = 0, where E^ = 
~ fr = 1230 — 502 MeV is the position of the A(1230) resonance, and a refers to the P33 
partial wave. This illustrates the fact that our quasi-particle propagator is the free NNA 
propagator for which the A has a width, and this width has the correct energy dependence 
as dictated by unitarity and the experimental phase shifts in this channel. Thus if we 
include the contribution to the three-body force from the P33 ttN channel only, we have 
effectively included, to lowest order, the contribution from the three-body force resulting 
from the A^A^ — A^A coupling. However, our result differs from the standard definition of 



this contribution to the three-body force |^5[ in that we have employed the nN data rather 
than the A^A^ data to fix the parameters of this force, and our A is a proper resonance 
in the vrA^ system and not a real particle as it is often considered. The other difference 
between this approach and that used in the A^A^ — A^A coupled channel approach is in 
the choice of the pion propagator in the "transition potential" . In the A^A coupled channel 
approach, the pion propagator is taken to be (/c^+m^)~^, which corresponds to the Feynman 
propagator with the nucleons on-mass- shell. This propagator has no energy dependence. On 
the other hand we have chosen the standard non-relativistic four particle Green's function 
with a relativistic expression for the pion kinetic energy. This is consistent with four-particle 
unitarity and is equivalent to a time ordered propagator, and to that extent we have only 
one time order for our pion exchange. Finally, we should note that the A^A^ — A^A coupled 
channel approach treats the contribution of the A to all orders and as a result includes the 
dispersive contribution which to a large extent cancels the contribution of the three-body 
force. 



The recent three- nucleon force results, reported by Pefia et al. |23], treat the A as a 
ttN resonance. To that extent, the vrA^ amplitude in the P33 channel is similar to that 
presented here, in that the mass and width of the A have energy dependence as dictated 
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by nN scattering data. Pefia et al. have the additional advantage that they not only have 
included the A contribution to the three-body force to all orders, but have also included 
the dispersive contribution. However, by restricting their Hilbert space to N and A and 
turning off all interaction in the pionic part of the Hilbert space, they have not included the 
contribution of the vrA^ amplitude in other than the P33 partial wave. They find that the 
contribution of the non-resonant vrA^ amplitude is very small as a result of the fact that the 
pionic component of the three-nucleon wave function is small. 

We now try to relate our three-body force with that used in the TM approach. In the 
TM approach, the final vrA^ amplitude is written in terms of the pion momenta, where as 
we use the ttN relative momenta Q'^ and Q3 [Q. To achieve their result, we recall that the 
partial wave expansion of the vrA^ amplitude is given by 

{Cl\t{E)m = Y.PM,^)UQ,Q':E) , (4.17) 

a 

where a = {£, j, T) are the quantum numbers corresponding to the orbital and total angular 
momentum and isospin of the tcN system, i.e. a labels the different partial waves. The 
partial wave projection operator Pq,(Q, Q') is defined as 

PaiQM') = PrT^i^iy^^nd {yirM') , (4.18) 

m 

with Prp the projection operator for a given isospin channel, and (Q|3^£jm) the eigenstates 
of the orbital and total angular momentum of the ttN system. If we now write the angular 
momentum projection operator in terms of Q ■ Q' and Q x Q' ■ cr and the isospin projection 
operator in terms of the Pauli isospin operator, we can write our vrA^ amplitude, assuming 
the pion scatters of nucleon 1, in the form presented by the Tucson-Melbourne formulation 
as 



F^'iQ, Q') = (t2 ■ T3) [a + bQ-Q' + ciQ'' + Q")J + rf(ri • {r, x T2))(o-i ■ (Q x Q')) . 

(4.19) 

For the factors a, b, c and d to be constant, as required by the Tucson- Melbourne definition 
of the three-body force, we have to make the following approximations: (i) Since the TM 
potential is derived from the off-mass-shell vrA^ amplitude T(i/, t, g, g'), where t is the Man- 
delstam variable and q {q') is the four momentum of the initial (final) pion, by expanding 
the amplitude about u = 0, we need to determine the corresponding approximation for our 
off-energy-shell nN amplitude. For the off-mass-shell amplitude the nucleon pole {s = m\i) 
traces a curve in the v — t plane that crosses the z/-axes close to v = Q ||5l|]. Since we 



find it difficult to directly relate the variables that in the off-mass- shell amplitude with the 
corresponding variables in off-energy-shell amplitude, we have chosen the position of the 
nucleon pole, i.e. E = rrijsi or s = m|r, to be the closest approximation to = 0. In this way 
we approximate the energy in our amphtude to be the nucleon mass, i.e., t{E2) — >■ r(mAr). 
(ii) The separable potential form factors QaiQ) have to be expanded in a power series in the 
momentum Q, keeping those terms such that the final amplitude does not have any powers 
of the momentum higher than the momentum squared. With these approximations, the 
Thomas separable potential, in conjunction with the potential PJ in the Pu channel, can 
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be written in the above form with the constants a, b, c and d given in Table |nT| . Included in 
the table are also the corresponding parameters from the Tucson-Melbourne potential |15 . 

In next section, we will present our numerical result and show: (i) How the contribution to 
the three-body force from each vriV partial wave depends on the energy in the nN amplitude, 
(ii) How the three-nucleon force is sensitive to the choice of the ttNN form factor which is 
determined as the residue of the Pn pole term. Though we have not included all the 
diagrams that would contribute to the three-body force, we have included the most important 
contribution, which will allow us to examine these effects. 



V. NUMERICAL RESULTS 

Having defined our three-nucleon force in terms of the ttN potential whose parameters 
have been adjusted to fit the vrA^ data, and in particular, the phase shifts up to the threshold 
for pion production, we turn our attention in this section to the calculation of the contribu- 
tion of this force to the binding energy of the three-nucleon system. As a first calculation 
with a three-body potential that includes the energy dependence of the vrA^ amplitude, we 
have chosen to use the first order perturbation theory to calculate this three-body force 
contribution. Therefore, we can write our three-body force as 

! ,j = l 

= E E w^.(p„q,;pU:) , (5.1) 

1 a 

where j in the sum refers to the nucleon that emits the pion, and i refers to the nucleon that 
absorbs the pion. Here, a refers to the tcN partial wave used to calculate the three-body 
force. Since the three-nucleon wave function used is a solution of the Faddeev equations, we 
can write this wave function as the sum of three components, i.e. 

1^) = l^i) + 1^2) + 1^3) 

= [1 + (123) + (132)] , (5.2) 

where we have written the total wave function in terms of the elements of the permutation 
operators. Making use of the properties of the permutation group, we can write the total 
contribution to the binding energy from this three-body force as 

AE(^) =6(^1 Vr(l, 3)1^) . (5.3) 

The three-body potential W{1, 3) can now be partial wave expanded in terms of the partial 
wave potential given in Eq. ( [4.13| ) as 

(Piqi|iy(l,3)|p3q3) =E E {pAi\Pi)Wg^^{puqi,P3,q3;E){f3,\p,q,) , (5.4) 

where (3i defines the three-body partial wave quantum number in which nucleon i is the 
spectator. In Eq. ( ^.4|) , a refers to the ttN partial wave that contributes to the three-body 



17 



potential. In a similar manner we have to expand the three-nucleon wave function in terms 
of the angular momentum and isospin bases (piqilA). This is given by 



/3i 



N2 

E (piqilA) 

/3i=l 



iVi 

i=2,3 I3j=l 



(5.5) 



In Eq. (|5.5| ) the sum over A^^i extends to the number of three-body channels included in the 
solution of the Faddeev equations which in turn is determined by the number of two-body 
A^A^ channels included. On the other hand, the sum over N2 is an infinity sum which we 
have truncated for practical calculations. For N2 > Ni the first term {piqilip/s-^) contributed 
only to the first A^^i terms in the N2 sum. This restriction is the result of truncating the 
number of A^A^ partial waves included in the solution of the Faddeev equations. We will 
test the convergence of our final results to both the sum over A^^i and A^2- Making use of the 
partial wave expansion for the three-body potential and the three-nucleon wave function, we 
can write the total contribution of the three-body force to the binding energy of the triton 
as 

00 CXD 00 CX) 

AE(3) = 6 J dPiPlJ dpspl J dq,ql j dq.ql 

" P-i-P'i 0000 
x(^|pigi;/?i)l^^%3(pi,gi,P3,g3;^) {h^qm^ , (5.6) 

where E is the energy of the three-nucleon system as determined by the solution of the 
Faddeev equations for a given two-nucleon interaction, and the partial wave three-body 
potential, W'^^^g, as given in Eq. ( [4.13[ ). 



In the present investigation we have chosen to use the Paris nucleon-nucleon poten- 



tial ||52|. The Paris potential is energy independent while our three-nucleon force has been 
derived to be energy-dependent. To that extent, our two- and three-body potentials are 
not consistent in that they are not derived from the same Lagrangian. However, since we 
are examining the energy dependence of the three-body force for each vrA^ partial wave, 
we hope that the present perturbative calculation may allow us to gain some insight into 
this problem. To simplify the construction of the three-nucleon wave function needed to 



evaluate the integrals in Eq. (|5.6| ), we have chosen the separable representation of the Paris 
potential ||^. This representation has been tested for the three-nucleon observables with 
considerable success ||53| . In Table 0, we give the rank of the separable expansion we have 



chosen. This choice was dictated by the requirement that we should reproduce the binding 
energy of the triton and the different percentages of S—, S'—, and D— state probability for 
the triton. However, before we compare our results with the coordinate space calculation for 
the Paris potential, we present in Table |V] the convergence of these quantities as we increase 
the number of three-body channels, A^i, in the solution of the Faddeev equations. From this 
table we can see that as far as the three-nucleon wave function is concerned, the 18 channel 
Faddeev equations give good convergence for all quantities. These 18 channels in the Fad- 
deev equations correspond to the truncation of the A^A^ interaction to include all two-body 
partial waves with total angular momentum less than or equal to two, including the coupled 
F2 NN channels. To justify the use of the separable expansion to the Paris potential. 
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PEST p2|], we compare, in Table 0, our results for the 18 channel calculation with the 
corresponding results based on the coordinate space solution of the Faddeev equations for 
the Paris potential. Taking the difference between the two coordinate space calculations as a 
measure of the numerical uncertainty in the solution of the Faddeev equations, we have good 
agreement with previous results for the Paris potential. This suggests that the three-nucleon 
wave function resulting from the PEST approximation is comparable to that resulting from 
a solution of the Faddeev equations in coordinate space for the exact Paris potential. 

Having established the fact that the three-nucleon wave function generated by the sep- 
arable expansion to the Paris potential is of comparable quality to that resulting from the 
solution of the Faddeev equations in coordinate space, we turn our attention to the conver- 
gence of the contribution of the three-body force to the number of three-body channels in 
the solution of the Faddeev equations, A^i, and the number of three-body channels included 
in the partial wave expansion of the wave function, N2- For this study we make use of the 



Pii potential PJ of McLeod and Afnan In Table |Vli| we present the contribution to 
the binding energy from the more important vrA^ partial waves for 5, 10, and 18 channel 
Faddeev calculations. In all cases we have taken 18 partial waves for the expansion of the 
wave function. All energies in Table [VI I| are in keV. We note at this stage that although 
the contribution to the binding energy is small, the 18 channel Faddeev calculation has 
converged, while the 5 channel calculation gives an incorrect result. With the 18 channel 
Faddeev calculation, we tested the convergence of our result to the number of terms in the 
partial wave expansion of the three-body wave function, A^2- From the results in Table [VIII 
we may conclude that N2 = 18 is sufficient to give us a 1 keV accuracy for the contribution 
from a given vrA^ partial wave. If the need arises we might have to resort to more terms in 
the partial wave expansion of the wave function. 

In Table |X| we present the contribution of the three-body force to the binding energy 
of the triton from the different ttN partial waves for two different Pu potentials. Here we 
have taken Ni = N2 = 18, with all energies given in keV. The most surprising result of 
our calculations is the overall small contribution of the three-body force. From the results 
in Tables VII and VIII , it is clear that the inclusion of more three-body partial waves in 
the wave function expansion and the solution of the Faddeev equation will not change the 
results substantially. Before we address the origin of this small contribution from our three- 
nucleon force it is interesting to note: (i) The comparable contribution from the S— and 
P— waves irN partial wave, and in particular the large contribution of the 6*31 compared to 
the 5*11. This suggests that we need to include both S— and P— wave nN amplitudes into 
the calculation. Furthermore, models based on the dominance of the A(1230) resonance 
might not be valid since they neglect the contribution from the S31 and Pu partial waves. 
In fact, in the present formulation, the contribution of the P33 partial wave amplitude is for 
ttN energies below the nucleon pole, some 300 MeV below the A resonance, (ii) There is a 
cancellation between the S— and P— wave nN contributions requiring a consistent treatment 
of both sets of partial wave amplitudes, (iii) The P— wave contribution comes equally from 
the non-pole part of the Pu, and the P33 partial wave amplitudes. This is despite the 
fact that the overall Pu phase shifts are small when compared with the P33 phase shifts. 
However, if we recall that it is the non-pole part of the amplitude that contributes to the 
three-body force, and this non-pole part, on its own, has phase shifts that are comparable 
to those in the P33 channel , then the results reported in Table [I^ are not surprising. 
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Finally, if we compare the results for the two potentials, we find that the potential Ml 
gives a larger contribution to the binding energy than the potential PJ. To understand this 



difference, we compare the dressed form factors for these two potentials in Fig. |10[ Here 
we observe that the potential Ml has a harder form factor than the potential PJ, i.e., the 
dressed form factor for potential Ml is greater than the corresponding form factor for the 
potential PJ for large k. This is consistent with results of the fact that the three-body force 
contribution to the binding energy increases as the form factor gets harder. We will come 
back to this point later in our discussion when we consider the role of the nNN form factor 
in the contribution of the three-nucleon force to the binding energy of the triton. 

We now turn to the question of why the contribution of this three-body force is small. 
From Table |T|, we may expect the maximum difference between 

our prediction and the TM result for the three-body force contribution to be at most 
an order of magnitude, but not three orders of magnitude. Since the unique feature of the 
present calculation is the inclusion of the energy dependence in both the irNN form factor 
and the nN amplitude, we will commence by turning this energy dependence off. We will 
also concentrate on those vrA^ partial waves that give a substantial contribution to the three- 
body force. We will restrict our results to Ni = N2 = 18. As a first approximation, denoted 
(i) in Table we fix the energy in the nNN form factor to be the position of the nucleon 
pole, i.e. 

f^ik;E)^f^ik;m,,) . (5.7) 

Although this approximation changes our final result by increasing the contribution of the 
three-body force to the binding energy, the magnitude of the increase is not substantial be- 
cause the cancellation between the repulsive 5*31 and the attractive Pu and P33 contributions 
is still present. In particular, we note that both the attractive and repulsive contributions 
have increased in magnitude. We next take the energy dependence in the irN amplitude to 
be the position of the nucleon pole (approximation (ii)), i.e. 

Ta{E) ^ T^{mN) . (5.8) 

In this case, we have increased the total contribution of the three-body force to the binding 
energy by an order of magnitude as compared to the exact result. This substantial increase 
in binding is mainly due to the fact that the 5*31 contribution is reduced in magnitude while 
the Pu and P33 contributions have increased, thus reducing the cancellation between the 
attraction and repulsion when compared with the exact calculation. To understand why 
the contribution of the S^i partial wave is suppressed as the energy in is increased (i.e. 
brought closer to the ttN threshold, see Fig. |^), we recall from Eqs. ( p.3|) and (|3.4| ) that 
A = +1 for repulsive potentials such as the 5*31, while A = —1 for attractive potentials such 
as the P33. On the other hand, (fi'a|G'^Ar(-E)|5'a) is negative for E < (m^ + rriTv) and increases 
in value as we approach the nN threshold from below. Thus for attractive potentials there 
is a cancellation in the denominator of Ta{E) giving rise to an increase in the value of 
Ta{E) as E approaches the threshold. On the other hand, for repulsive potentials, e.g. the 
5*31 channel, the value of Ta{E) decreases as we approach the elastic threshold from below, 
resulting in a suppression of the repulsive contribution. To demonstrate the validity of this 
argument, we have proceeded to change the energy in the vrA^ amplitude, to be the threshold 
for ttN scattering, i.e.. 
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TaiE) -> Ta{mN + m^) . (5.9) 

The result of this approximation is labeled (iii) in Table |^. Here, we observe that the 
contribution of the S^i is further reduced, while the Pu and P33 contributions are increased in 
magnitude, giving even more attraction. We now fix the energy of both the uNN form factor 
and the vrA^ scattering amplitude to the nucleon pole, i.e. E = rriN, (i) & (ii) in Table 
This gives even more binding than fixing the energy in either the ttN scattering amplitude 
or the TiNN form factor. Finally, we can increase binding further by fixing the energy of 
the nNN form factor at the nucleon pole, while the energy in the vrA^ amplitude is taken to 
be the threshold energy, (i) & (iii) in Table From the above analysis we may conclude 
that it is the energy dependence in both the nNN form factor and the nN amplitude that 
has substantially reduced the contribution of the three-nucleon force to the binding energy 
of the triton, and this reduction is a result of the cancellation between the repulsive 5*31 
contribution and attractive Pu and P33 contributions. However, the approximation of fixing 
the energy in the vrA^ amplitude and ttNN form factor does not give us a sufficiently large 
contribution which is comparable with the result for the TM potential. 

Finally, to fully understand the origin of this small contribution to the binding energy 
from the three-nucleon force, we turn our attention to the form factors used in the separable 
potential and the form factor in the tcNN vertex. It is common practice to take the ttNN 
form factor in A^A^ scattering to be either of a dipole or monopole form. Therefore, as a first 
step in examining the sensitivity of our final three-body force contribution to the binding 
energy, we replace the nNN form factor f^{k; E) by a monopole, i.e. 



k 

where the monopole form factor -Fo(^) is given by 



kFoik) , (5.10) 

fc=0 



foW = ^ (5.11) 



with the cutoff mass A, varied. In Table ^ we compare exact results for the Pu potential 



PJ with the result for the approximations in Eqs. ( [5.7| ) and ( |5.8| ), referred to as (i) & (ii 



and the approximations in Eqs. ( |5.8|) and (|5.10|) with A = 400 and 800 MeV. The results 
in lines 2, 3 and 4 of Table ^ have the energy in the vrA^ amplitude fixed to be the energy 
at the nucleon pole, i.e. E^^j^ = mjy. Here we observe that the final total contribution 
of the three-body force increases by an order of magnitude when the energy in both the 
nNN form factor and nN amplitude is fixed at the nucleon pole. There is a further order 
of magnitude increase when the irNN form factor is replaced by a monopole form factor 
with a cutoff mass of 400 and then 800 MeV. In fact, half of this increase is achieved when 
the cutoff mass is increased from 400 to 800 MeV. This establishes the sensitivity of our 
result to the choice of tcNN form factor. In this case it is interesting to observe that the 
5*31 and Pn contributions to the three-nucleon force almost completely cancel, leaving the 
P33 contribution to be approximately the total contribution. Thus in this approximation, 
the contribution to the three-nucleon force is predominantly due to the channel in which 
the A(1230) dominates the scattering amplitude. This result should be compared with the 



exact results, line 1 of Table where the S31 contribution cancels the sum of the Pu and 
P33 contributions. 
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In the spirit of the TM approach, the ttN amphtude is expanded to lowest order in ^ 
and a monopole ttNN form factor is introduced. This approximation may be implemented 
by the following replacement: 



TaiE) 

9aik) 



k 

Ta{niN) 

9aik) 



kFoik) 



k=0 



(5.12) 



k^ 



k'Foik) 



k=0 



where i is the angular momentum in channel a and the monopole form factor Fo{k) is 
defined in Eq. ( |5.11|) . However in our formulation, this replacement will destroy the fit to 



the experimental vrA^ phase shifts. In Table pCIl] , we compare our exact result for the Pn 
potential PJ |^^, with the results of the approximation in Eq. (|5.12| ) with A = 400 and 
800 MeV. The approximation in Eq. ( ^.121) gives rise to an increase in the contribution 
of the P-waves substantially, while the S-wave contribution remains relatively unchanged, 
and therefore negligible. In fact we can now adjust the cutoff mass A to get the difference 
between the experimental binding energy and the calculated three-nucleon result for any of 
the two-nucleon interactions. 

To understand this large change in the magnitude of the total contribution of the three- 
body force when we introduce the form factor Folk), we compare in Fig. the monopole 
form factor with A = 400 and 800 MeV and the dressed nNN form factor for the potentials 
PJ and Ml. Here we observe that the form factor -Fo(^) with A = 800 MeV is almost a factor 
of 3 larger than the dressed irNN form factor at ~ 3 fm~^. Furthermore, this form factor 
comes raised to the power of two in the case when only the ttNN form factor is replaced 
by a monopole form factor, or a power of four when the approximations in Eq. ( |5.12|) are 
implemented, i.e. we have a power of two from the tcNN vertices and another power of two 
from the form factor of separable ttN amplitudes. Thus the difference between the result of 
including f^{k, tjin) and -Fo(^) should be roughly one order of magnitude in Tables |X|, and 



two orders of magnitude in Table pCII| . In the latter case we assume the separable potential 
form factors have similar ranges to the dressed ttNN form factor. This explains the drastic 
change in the contribution of the three-body force to the binding energy when we introduced 
the monopole form factor -Fo(^) into our calculation. 

Finally, to get the closest approximation to the TM three-nucleon force, we have modified 
the propagator for the exchanged pion by replacing our propagator by the corresponding 
Feynman propagator, i.e.. 



Qo — 



2tOr 



1 



Q"^ - ml 



2uJr 



(5.13) 



where uq = yQ'^ + ml.. To understand the difference between these two propagators, 
we recall that the Feynman propagator has both a positive energy and negative energy 
component, since 



1 



1 



mt 



2uor 



1 



1 



(5.14) 
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and our choice of propagator, i.e. (Qo — ^q)~^7 corresponds to taking the positive energy 
component of the Feynman propagator. However in any energy independent approximation, 
the Feynman propagator reduces to —I/uJq. The substitution in Eq. ( |5.13| ) gives us the last 
hne of Table PCII] , and results in the reduction of the contribution of the three-body force 
so that a monopole form factor with a cutoff mass of 800 MeV will give a three-body force 
contribution of about 0.6 MeV, which is consistent with the results reported in the literature 
for the TM three-nucleon force. 

Thus, to get a substantial contribution from the three-body force we have had to make 
two approximations, (i) We have dropped the energy dependence of the ttN amplitude and 
the ttNN form factor, (ii) We have modified the off-shell behavior of the ttN amplitude 
by introducing the same monopole form factor in all partial waves, at the sacrifice of the 
fit to the experimental data, in order to get a substantial increase to the three-body force 
contribution. Although the first approximation is not justified, the second could be accepted 
on the ground that we have chosen the wrong off-shell behavior. In particular, we should 
consider changes the cutoff mass in the bare nNN form factor fo{k), and make use of the 
Goldberger-Treiman relation to constrain the dressed ttNN form factor f^{k] E). This 
point is presently under further investigation. 



VI. CONCLUSION 

The main motivation for this investigation was to establish the relative contribution of 
the different partial waves of the ttN amplitude in determining the vr — vr three-nucleon 
force and its contribution to the binding energy of the triton. To achieve this, we made 
use of a separable potential formulation of the irN scattering, taking into consideration 
that such a parameterization of the off-shell tiN amplitude has been used in ird scattering 
and in the derivation of the pion-nucleus optical potential with considerable success. In 
particular, we maintained the energy dependence of the nN amplitude since that could 
effect the overall contribution to the three-body force. More importantly, the cancellation 
between the contribution of the different partial waves could be sensitive to the inclusion 
of this energy dependence. Since the energy domain, important to the determination of the 
three-nucleon force, is in the unphysical region and below the position of the nucleon pole in 
the energy plane, we chose to fit the scattering data closest to this region. In particular we 
chose to fit the scattering lengths, the phase shifts up to the pion production threshold, and 
the position and the residue of the vrA^ amplitude at the nucleon pole. This parameterization 
allowed us to determine the ttNN form factor and that part of the vrA^ amplitude that gives 
rise to the three-body force. 

From the numerical results of our calculation we can conclude that: (i) The energy 
dependence of both the nNN form factor and vrA^ amplitude gives a suppression of the 
contribution of Fig. ^ to the three-nucleon force. Given the fact that this diagram has 
always been considered to give the main contribution to this force, we can conclude that the 
three-nucleon force for this ttN parameterization is small, and will not change substantially 
if we include this three-body force in an exact calculation rather than in the perturbative 
approach used in the present investigation, (ii) The inclusion of the energy dependence gives 
rise to a substantial cancellation between the contribution from the repulsive 5*31 partial wave 
and the attractive Pu and P33 partial waves, (iii) The contribution from the P33 partial 
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wave is not as dominant as we would expect. In fact the attraction comes equally from 
the Pii and P33, while the repulsion comes from the S31. This raises a question about 
the validity of including the three-body force in terms of the A^A^ — A^A coupled channel 
approach while neglecting the S"— wave and the Pn tvN contributions, (iv) The choice of the 
ttNN form factor, to be determined by the residue of the nN amplitude at the nucleon pole, 
is the other main reason for the reduction in the overall magnitude of the three-nucleon force 
contribution to the binding energy. The question of the possibility of choosing the cutoff 
mass in the bare form factor fo{k) to be consistent with the ttN data, while maintaining a 
substantial three-body force, will need further investigation. Furthermore, the Goldberger- 
Treiman relation p!3| , p!^ should be used to constrain the momentum dependence of the 
dressed irNN form factor f^{k; E). 

To further substantiate the above conclusions, we may need to examine a number of 
questions: (i) How sensitive are our results to the choice of the nN interaction? In particular, 
would we get a small contribution from the three-nucleon force, and specifically the diagram 
in Fig. |I|, if we commenced with a chiral Lagrangian such as that used by Pearce and 
Jennings [0? (ii) Would our final results be substantially different if this three-body force 
were to be included to all orders? Examination of the perturbation series for the Paris 
potential with the TM three-body force [^] suggests that higher order contributions are 
not negligible, but the magnitude of the overall contribution does not change by more than 
a factor of two. This further suggests that an exact calculation will not effect our final 
conclusion, (iii) How important are the dispersive effects? The latest results based on the 
A^A^ — A^A coupled channel approach suggests that there is a substantial cancellation 
between the three-body force contribution and the dispersive effects in the P33 channel. 
What happens to this cancellation when other vrA^ partial waves are included? If this 
cancellation is present for all vrA^ partial waves then it would lead to further reduction in 
the three-body force. 
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FIGURES 



FIG. 1. The contribution to the three-nucleon force. 

FIG. 2. The NN complex energy plane with the unitarity cut and the domain of integration 
when used in the Faddeev equation for the three-nucleon system. 

FIG. 3. The ttN complex energy plane with the unitarity cut, the N and A pole position, and 
the region of integration when calculating the three-body n — it force. 

FIG. 4. The dressing of nucleons with the restriction on the Hilbert space to include nucleons 
and up to one pion only. 

FIG. 5. The full dressing of nucleons on both sides of the ttA'^A'" vertex. 

FIG. 6. The lowest order terms resulting from the coupling of the NNN to the ttNNN channels 
as given in Eq.(2.17). (a) One pion exchange term, Vqpe- (b) Dispersion term, Vdisp- And (c) 
three body force, V^b- 

FIG. 7. The nucleon pole contribution to the ttN potential. 

FIG. 8. The contribution of the crossed diagram and the contact diagram to the ttN potential. 

FIG. 9. An illustration of the Jacobi momenta defined for the three-body force. 

FIG. 10. Comparison of the dressed ttNN form factors for the potentials PJ, and Ml with a 
monopole form factor having cutoff mass of A = 400, 800 MeV. 
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TABLES 



TABLE I. Parameters of the McLeod-Afnan Pn potentials. 

Pot. no rii Co ci C2 a j3i (32 rriQ 

fm~^ fm~^ fm^^ fm~^ 

PJ 2 i 43.5646 0.2907 412.97 3.8206 1.2688 5.181 5.1574 

A/1 1 2 1.0726 0.3433 2.54 2.7703 1.4422 2.1982 5.4314 



TABLE IL The scattering volume an, wave function renormalization Z2, and the irNN cou- 
pling constant at A; = for the Pu potentials. The fnNNiO) is to be compared with a value of 
0.2726 predicted by the Goldberger-Treiman relation of 3% change in fnNN{k) between k = ko 
and k = 0. 



Pot. 



On Z2 firNNiO) 

-3\ 



[m 



PJ 

Ml 



-0.0706 0.8059 0.2280 

-0.0721 0.7273 0.2209 



TABLE in. The parameters of the Tucson-Melbourne type irN amplitude as reported in Ref. 
1 15 1, and those extracted from the potential of Thomas and the Pn potential PJ |44|. 



Pot. 


a 


b 


c 


d 


Pefia & Coon 


1.03 


-2.62 


0.91 


-0.75 


Thomas &: PJ 


0.46 


-4.66 


0.67 


-2.47 



TABLE IV. The rank of the separable expansion to the Paris potential for the two-body NN 
partial waves included in the calculation of the wave function. 





^5*0 




'^Po,^Pi,^Pi 


^P2-^F2 


^D2,^D2 


Rank 


3 


4 


2 


3 


2 
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TABLE V. The convergence of the binding energy of the triton, the S— S'— and I?— state 
probabihties for 5, 10 and 18 channel Faddeev calculations. In calculating the wave function we 
have taken = 146. 





B.E. (MeV) 


P{S)% 


P{S')% 


P{D)% 


5 


7.2659 


89.882 


1.652 


8.401 


10 


7.0966 


90.274 


1.479 


8.182 


18 


7.3175 


90.111 


1.430 


8.393 



TABLE VI. Comparison of our 18 channel results for the separable expansion (PEST), with 
the results of coordinate space calculations of Los Alamos (LA) and Sendai (S) [Q]. 



Model 


B.E. (MeV) 


PiS)% 


P{S')% 


PiD)% 


PEST 


7.318 


90.111 


1.430 


8.393 


LA 


7.388 


90.130 


1.395 


8.409 


S 


7.56 


90.17 


1.32 


8.45 



TABLE VII. Convergence of the contribution of the three-nucleon force to the binding energy 
of the triton from different vrA^ partial waves and for different numbers of the three-body channels 
in the Faddeev equations. Here we have taken = 18. All energies are in keV. 



Ni 


5*11 


S31 


Pn 


P33 


5 


- 5.5 


31.8 


- 7.7 


- 15.2 


10 


- 4.8 


26.5 


- 7.9 


- 15.0 


18 


- 4.8 


26.4 


- 8.8 


- 16.0 



TABLE VIII. Convergence of the contribution of the three-nucleon force to the binding energy 
from different vrA^ partial waves as a function of A'2, the number of three-body channels in the 
partial wave expansion of the wave function. The Faddeev equations are solved with 18 channels. 
All energies are in keV. 



N2 


5*11 


'S'31 


Pn 


^33 


10 


- 5.2 


37.1 


- 7.2 


- 8.5 


18 


- 4.8 


26.4 


- 8.8 


- 16.0 


26 


- 5.3 


26.2 


- 9.5 


- 16.1 


34 


- 5.3 


25.7 


- 9.5 


- 16.9 
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TABLE IX. The contribution to the three-body force from different vrA'^ partial wave amph- 



tudes. The results are for two different choices of the Pu potential. All energies are in keV. 



Pot. 




Sn 


^31 


Pii 


P31 


Pl3 


-P33 


Total 


PJ 




-4.8 


26.4 


-8.8 


-0.0 


4.5 


-16.0 


-2.3 


Ml 




-5.0 


28.9 


-15.3 


-2.1 


6.2 


-22.1 


-9.4 


TABLE X. The effect of removing the energ 


nr r\ (^tt^ tin r\ t^n fPi 


in the ttNN form factor 


(i), and 


the ttN amplitudes, (ii] 


and (iii). 


The results in 


llilO IdUit/ dlt/ 


for the Pii 


potential PJ, and all 


energies arc in keV. 














Approx. 


Su 


Ssi 


Pn 


^ 61 


^^13 


-F33 


Total 


exact 


-4.8 


26.4 


-8.8 


-3.6 


4.5 


-16.0 


-2.3 


( 





-5.9 


31.7 


-13.0 


-3.0 


5.7 


-22.3 


-6.8 


(ii) 


-5.0 


23.5 


-13.7 


-3.6 


3.9 


-28.2 


-23.1 


(iii) 


-5.7 


19.5 


-19.6 


-3.5 


3.6 


-49.8 


-55.5 




. (ii) 


-6.1 


28.1 


-20.1 


-3.1 


5.0 


-39.8 


-36.0 


(i) & (iii) 


-7.0 


23.4 


-28.7 


-3.0 


4.7 


-70.3 


-80.9 



TABLE XI. The effect of changing the ttNN form factor on the contribution of the 
three-nucleon force to the binding energy of the triton. The energy in the ttN amplitude is fixed 
at m-N- The comparison is between monopole form factor with a cutoff mass A = 400 or 800 MeV 
and f^{k, niN)- Also included are the exact results which have the full energy dependence of both 
the ttNN form factor and nN amplitude. The total includes the contributions from all S- and 
P-wave ttN amplitudes. All energies are in keV. 



A 


^11 


^31 


Pii 


-P33 


Total 


exact 


-4.8 


26.4 


-8.8 


-16.0 


-2.3 


(i) & (ii) 


-6.1 


28.1 


-20.1 


-39.8 


-36.0 


400 


-12.5 


57.6 


-40.7 


-81.2 


-68.7 


800 


-29.3 


147.5 


-204.1 


-364.0 


-395.1 
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TABLE XII. The effect of changing the form factors on the contribution of the three-nucleon 
force to the binding energy of the triton. Here, both the ttNN form factors/^(fc; ii^) and the 
separable potential form factors gaik) are replaced by a monopole with a cutoff mass A = 400 or 
800 MeV. The last line in this table corresponds to taking the 'Feynman' propagator for the pion. 
The total includes the contributions from all S- and P-wave ttN amplitudes. All energies are in 



keV. 



A 


5ii 


•S'si 


Pu 


-P33 


Total 


exact 


-4.8 


26.4 


-8.8 


-16.0 


-2.3 


400 


-15.0 


18.8 


-58.8 


-211.5 


-231.1 


800 


-82.7 


113.2 


-1146.7 


-3172.5 


-3897.2 


800* 


-36.0 


47.4 


-193.4 


-542.5 


-611.6 
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